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ABSTRACT 


The  research  described  in  this  thesis  is  a  continuation  of  work  started  by  the 
Applied  Research  Laboratories  of  the  University  of  Texas  at  Austin  into  the  analysis  of 
biosonar  signals.  Experiments  conducted  in  1997  on  two  species  of  small  toothed 
whales,  found  these  species  to  emit  significant  high  frequency  signal  components, 
extending  to  as  high  as  400  to  500  kHz. 

To  assess  the  importance  of  these  high  frequencies  in  dolphin  echolocation  and 
target  identification,  experiments  were  performed  in  which  an  acoustic  filter,  used  to 
suppress  the  high  frequencies,  was  placed  between  a  dolphin  and  a  target.  Insertion  Loss 
and  Reflection  Loss  measurements  performed  on  14"  thick  and  14"  thick  Sound  Absorbing 
Filters  (SOAB)  demonstrated  their  effectiveness  at  absorbing  high  frequencies  above  150 
kHz,  with  little  reflectivity. 

The  results  from  one  echolocation  experiment,  with  one  dolphin,  showed  the 
animal's  ability  to  classify  targets  was  essentially  unaffected  by  the  insertion  of  the  filters. 
Analysis  of  the  dolphin's  echolocation  signals  showed  the  animal  definitely  compensating 
for  the  filters,  by  increasing  its  sound  energy  output,  especially  at  frequencies  above  100 
kHz.  It  is  anticipated  that  this  initial  experiment  will  lead  to  future  research  in 
explaining  the  existence  of  these  high  frequency  echolocation  components. 
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I.  INTRODUCTION 


The  work  described  in  this  thesis  is  a  continuing  research  effort  started  by  the 
Applied  Research  Laboratories  of  the  University  of  Texas  at  Austin  (ARL-UT),  and 
sponsored  by  the  Office  of  Naval  Research,  into  the  analysis  of  biosonar  signals. 
Researchers  at  ARL-UT  have  identified  significant  high-frequency  components  in 
dolphin  echolocation  signals  [Ref.  1].  The  primary  focus  of  the  present  research  is  to 
further  examine  these  signals,  and  to  conduct  blind  dolphin  biosonar  target  detection 
experiments  with  and  with-out  high-frequency-absorbing  screens.  These  experiments 
will  assess  the  importance  of  these  high  frequencies  in  dolphin  echolocation  and 
identification. 

A.  BACKGROUND 

Over  the  past  35  years,  most  acoustic  experiments  performed  with  dolphins  have 
utilized  hydrophones  with  a  typical  receiving  sensitivity  curve  extending  to  only  about 
130  kHz.  As  a  result,  prior  measurements  on  the  dolphin  echolocation  signal  typically 
show  the  main  frequency  components  peaking  at  around  100  kHz.  The  prior  work  shows 
that  the  signal  decays  at  frequencies  approaching  the  hydrophone  maximum  receive 
sensitivity,  typically  around  130  to  150  kHz.  [Ref.  2] 

In  a  project  summary  report  to  the  Office  of  Naval  Research  Biosonar  Research 
Program  entitled  “Role  of  Nonlinear  Acoustics  in  Biosonar,”  Professor  Thomas  Muir  and 
his  graduate  student,  Ms.  Diane  Blackwood  of  Texas  A&M,  summarized  the  results  of 
experiments  conducted  in  May  1997.  These  experiments  were  conducted  at  the  Naval 
Command  and  Control  and  Ocean  Surveillance  Center  (NRaD),  in  San  Diego,  California, 
and  involved  the  bottlenose  dolphin  ( tursiops  truncatus )  as  well  as  the  beluga  whale 
(delphinapterus  leucas).  Their  experiments  proved  that  these  species  do  emit  higher 
frequency  echolocation  signal  components,  extending  up  to  some  four  to  five  times  what 
previous  investigators  have  recorded. 
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Prior  research  on  the  bottlenose  dolphin  species  has  also  reported  the  high 
resolution  capability  of  marine  biosonars  to  classify  and  distinguish  between  small  man 
made  targets.  This  reported  level  of  resolution  could  not  be  achieved  with  man  made 
sonars  operating  in  the  frequency  range  around  100  kHz.  Cetacean  biosonar  performance 
reported  in  the  literature  has  sometimes  seemed  to  be  in  violation  of  a  law  in  physics 
called  the  “Uncertainty  Principle.”  This  principle  states  that  the  best  resolution  an  active 
sonar  beam  can  achieve  equals  the  sonar  pulse  duration  times  the  medium  sound  speed, 
divided  by  two,  (cx/2)  [Ref.  3].  Using  this  relation,  it  is  predicted  that  a  dolphin 
projecting  an  individual  click  of  60  ps,  at  a  hollow  steel  cylinder  target,  and  using  1500 
m/s  for  the  sound  speed  of  seawater,  can  achieve  a  best  resolution  of  4.5  cm.  Previously 
reported  biosonar  experiments  on  dolphins  have  shown  they  are  capable  of  achieving 
much  higher  resolution.  The  capability  of  a  bottlenose  dolphin  to  discriminate 
differences  in  the  wall  thickness  of  hollow  steel  cylinders  was  studied  by  Titov  [Ref.  4]. 
The  animal  was  able  to  react  to  a  wall  thickness  difference  of  0.2  mm  at  a  75%  correct 
response  level.  The  existence  of  high  frequency  signal  components  (greater  than  100- 
150  kHz)  may  help  to  explain  why  a  dolphin  can  achieve  a  much  higher  resolution  than 
permitted  by  this  law  in  physics.  [Ref.  1] 

The  dolphin’s  ability  to  recognize  and  classify  targets  buried  in  the  sediments,  in 
reverberation  limited  environments,  is  better  than  any  man-made  mine-hunting  sonar 
system.  In  fact,  marine  mammals,  although  cumbersome,  and  expensive,  are  currently 
the  only  means  the  Navy  has  for  detecting  buried  mines  [Ref.  5].  Therefore,  a  brief 
description  of  the  current  U.S.  Navy  Marine  Mammal  Program  is  given  in  section  C. 

B.  RESEARCH  MOTIVATION  AND  OBJECTIVES 

This  thesis  describes  the  results  of  an  investigation  into  the  effect  of  the  insertion 
of  a  filter  which  suppresses  the  high  frequencies  on  the  dolphin’s  ability  to  classify 
targets.  This  thesis  research  has  significance  to  both  the  military  and  commercial 
interests.  For  the  military,  the  understanding  of  biosonar  mechanisms  and  signal 
production  can  be  utilized  in  improving  the  resolution  of  U.S.  Navy  sonar  systems,  which 
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have  a  much  higher  area  search  rate  than  marine  mammals.  An  understanding  of  the 
mechanisms  that  enable  dolphins  to  detect  buried  mine-like  objects  could  lead  to 
considerable  improvements  in  man-made  sonar  systems  for  buried  mine  detection  and 
classification.  This  knowledge  would  also  greatly  improve  the  biological  and  physical 
modeling  of  animal  acoustic  systems. 

C.  U.  S.  NAVY  MARINE  MAMMAL  PROGRAM 

The  Navy’s  Marine  Mammal  Program  incorporates  specially  trained  Atlantic  and 
Pacific  bottlenose  dolphins,  white  whales,  and  sea  lions  for  mine  detection  and 
neutralization,  swimmer  defense,  and  recovery  of  exercise  mines  and  torpedoes.  Taking 
advantage  of  years  of  evolution  that  have  produced  animals  well  suited  for  these  tasks, 
the  Navy  has  evolved  complex  and  sophisticated  training  techniques  that  enable  these 
animals  to  conduct  real-world  operations.  [Ref.  5] 

The  Marine  Mammal  Program  began  in  1960,  when  several  dolphins  were  used  in 
hydrodynamic  studies  addressing  underwater  torpedo  design.  In  1963,  the  Navy  began 
studying  the  animals’  deep  diving  and  echo-location  capabilities,  and  determined  that 
dolphins  could  work  untethered  in  the  open  ocean.  In  the  late  1960’s  the  Navy  developed 
a  dolphin  swimmer  detection  and  marking  system  under  the  code  name  Short  Time.  It 
deployed  to  Cam  Rahn  Bay  in  1970,  to  guard  an  ammunition  pier  that  had  been  the  target 
of  attacks  by  the  Vietcong.  Once  the  dolphins  were  on  scene,  the  raids  stopped.  In  1987, 
six  Pacific  bottlenose  dolphins  provided  underwater  surveillance  and  detection  capability 
to  support  bases  in  the  Persian  Gulf.  [Ref.  5] 

The  Navy’s  operational  Marine  Mammal  System  includes  four  to  eight  marine 
mammals  which  can  be  easily  deployed  on  very  short  notice  by  strategic  airlift  to  any  part 
of  the  world  and  can  be  worked  from  ships  in  forward  areas.  The  system  is  divided  into 
four  programs  utilized  by  the  fleet,  three  of  which  include  bottlenose  dolphins: 

•  Mk  4  Mod-0  -  Pacific  bottlenose  dolphins  detect  mines  and  attach 

neutralization  charges  on  the  mooring  cables  of  tethered  mines  moored  near 
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the  bottom.  The  Navy  is  expanding  this  system’s  capability  to  neutralize  all 
tethered  buoyant  mines. 

•  Mk  6  Mod-1  -  Dolphins  provide  defense  of  harbors,  anchorages,  and 
individual  ships  against  swimmers  and  divers.  The  Mk  6  participates 
regularly  in  fleet  exercises  and  real-world  base  security,  providing  a 
comprehensive  surface  and  subsurface  swimmer  detection. 

•  Mk  7  Mod-1  -  Dolphins  detect,  locate,  and  mark  or  neutralize  bottom  mines 
and  buried  mines.  This  animal  system  represents  the  only  operational  buried- 
mine  detection  and  neutralization  capability  in  the  world  today. 

The  Mk4  and  Mk7  Marine  Mammal  System  detachments  are  integral  operational 
elements  of  the  Navy’s  mine  countermeasures  forces  and  have  demonstrated  the 
capability  to  operate  for  extended  periods  from  ships  forward  deployed.  [Ref.  6] 

There  is  also  an  additional  system  under  development;  Experimental  8  Marine 
Mammal  System  will  employ  six  dolphins  for  exploration  and  reconnaissance  of  in¬ 
volume  moored  and  bottom  mine-like  contacts  in  the  Very  Shallow  Water  Zone  (10-40 
foot  depth).  The  Ex  8  dolphins  will  be  deployable  from  an  Amphibious  Task  Force  ship 
for  low- visibility,  minefield  exploration  and  reconnaissance  [Ref.  7]. 

The  dolphins  in  the  Marine  Mammal  Program  satisfy  critical  requirements  and 
real  world  operational  needs  that  today  cannot  be  met  as  effectively  or  efficiently  in  any 
other  way. 

D.  THESIS  OUTLINE 

The  second  chapter  provides  a  description  of  the  dolphin  echolocation  system  and 
characteristics  of  biosonar  signals  recorded  with  a  wide  band  hydrophone.  The  third 
chapter  describes  the  theory  and  laboratory  experiments  conducted  on  absorptive  acoustic 
screens.  The  fourth  chapter  explains  the  procedure  and  initial  results  using  an  acoustic 
filter  to  suppress  high  frequencies  and  its  effects  on  the  dolphin’s  ability  to  classify 
targets.  The  final  chapter  provides  concluding  remarks  and  recommendations  for 
continuing  research  efforts. 
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II.  DOLPHIN  ECHOLOCATION  SYSTEM 


This  chapter  will  present  a  brief  introduction  to  the  dolphin  biosonar  transmission 
system  and  the  characteristics  of  its  biosonar  signals.  In  addition,  samples  of  bottlenose 
dolphin  click  trains  will  be  analyzed  to  show  that  previously  ignored  and  undetected  high 
frequency  echolocation  signal  components  are  indeed  present. 

A.  DESCRIPTION  OF  THE  DOLPHIN  ECHOLOCATION  SYSTEM 

The  term  echolocation  refers  to  an  ability  the  dolphin  possesses  that  enables  it  to 
"see"  by  listening  for  echoes.  Figure  2.1  illustrates  the  echolocation  process.  The 
dolphin  echolocation  system  is  a  highly  specialized  sonar  that  enables  dolphins  to  explore 
their  environment  and  search  out  their  prey  in  a  watery  world  where  sight  is  often  limited 
by  dark,  murky  water  cluttered  with  debris.  How  a  dolphin  produces  and  receives  sound 
is  still  a  highly  controversial  subject.  Professor  Ridgway  proposed  a  predominant  theory 
that  the  nasal  plugs,  under  muscular  control,  produce  sound  in  the  form  of  acoustic 
transients  as  air  passes  between  the  plugs  and  the  nasal  walls  [Ref.  8].  The  frequency 
range  of  these  “sonar  clicks”  is  higher  than  that  of  the  sounds  used  for  communication, 
and  differs  between  species.  One  current  hypothesis  is  that  this  sound  is  projected  into 
the  water  in  a  narrow  beam  after  passing  through  a  fatty  melon  which  may  act  as  a  device 
to  couple  sounds  produced  deep  in  the  skull  into  the  water  [Ref.  9]. 


Figure  2.1  Dolphin  Echolocation  Mechanism  After  Ref.  [10]. 
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When  the  sound  strikes  an  object,  some  of  the  energy  of  the  sound  wave  is 
scattered  back  towards  the  dolphin.  It  has  been  postulated  that  sound  waves  in  the  water 
pass  into  the  head  of  the  animal  and  are  transmitted  to  the  ear  region  by  a  thin  bony  area 
in  the  panbone  within  the  dolphin’s  lower  jaw  [Ref.  9]. 

The  time  lapse  between  click  and  echo  could  enable  the  dolphin  to  evaluate  the 
distance  between  it  and  the  object,  as  is  the  case  in  torpedo  sonar,  for  example.  Professor 
Au  has  speculated  that  the  strength  of  the  signal  as  it  is  received  on  the  two  sides  of  the 
dolphin’s  head  may  enable  it  to  evaluate  direction  or  localize  sound  [Ref.  2].  By 
continuously  emitting  clicks  and  receiving  echoes  in  this  way,  the  dolphin  can  track  and 
find  objects. 

The  echolocation  system  of  the  dolphin  is  extremely  complex.  Using  only  its 
acoustic  senses,  a  dolphin  can  discriminate  between  practically  identical  objects,  which 
differ  by  ten  per  cent  or  less  in  volume  or  surface  area.  It  can  do  this  in  a  noisy 
environment,  can  whistle  and  echolocate  at  the  same  time,  and  can  echolocate  on  near 
and  distant  targets  simultaneously. 

B.  DOLPHIN  SONAR  RESEARCH  PRIOR  TO  1997 

Dolphins  are  capable  of  producing  extremely  short  duration,  broad  bandwidth, 
acoustic  signals,  which  are  utilized  for  echolocation.  The  ability  of  dolphins  to 
accurately  perceive  their  environment  and  to  perform  difficult  recognition  and 
discrimination  tasks  depends  on  the  characteristics  of  these  biosonar  signals  and  how  they 
are  emitted,  and  processed  upon  reception.  Signal  characteristics  and  projection  patterns 
have  been  recorded  and  studied  over  a  long  period  of  time  by  many  investigators,  but  the 
operational  mechanisms  of  dolphin  sonar  yet  remains  unanswered. 

A  typical  biosonar  signal  waveform  and  frequency  spectrum  of  a  bottlenose 
dolphin  recorded  in  a  tank  environment  by  Evans  in  1973  is  shown  in  Figure  2.2 
[Ref.  11].  The  peak  frequency  (frequency  of  maximum  energy)  in  this  example  was  52 
kHz.  Early  bottlenose  dolphin  signals  were  measured  in  tanks,  and  it  was  generally 
believed  that  peak  frequencies  occurred  in  the  vicinity  of  30  to  60  kHz. 
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In  1974,  Au  observed  significant  energy,  up  to  the  limit  of  his  detection  system, 
within  a  dolphin  click.  This  energy  extended  to  much  higher  frequencies  than  were 
previously  measured.  He  conducted  target  detection  experiments  in  Kaneohe  Bay,  Oahu, 
Hawaii,  which  involved  measuring  two  bottlenose  dolphin  echolocation  signals  in  open 
waters.  His  results  showed  that  the  signals  had  peak  frequencies  between  120  and  130 
kHz,  which  were  over  an  octave  higher  than  the  peak  frequencies  recorded  by  Evans. 

The  average  waveform  and  frequency  spectrum  of  a  biosonar  click  train  observed  by  Au 
is  shown  in  Figure  2.3.  [Ref.  2] 


UJ 


a. 


FREQUENCY  (kHz) 

Figure  2.2  Typical  Waveform  And  Frequency  Spectrum  of  Bottlenose  Dolphin  in 
a  Tank  From  Ref.  [1 1]. 


FREQUENCY  (  KHZ  ) 


Figure  2.3  Average  Waveform  and  Frequency  Spectrum  of  Bottlenose  Dolphin  in 
Open  Waters  From  Ref.  [2]. 


Mitson  published  evidence  of  high-frequency  acoustic  emissions  from  a  school  of 
white  beaked  dolphin  ( lagenorhynchus  albirostris )  in  the  North  Sea  in  1987.  While 
onboard  a  British  fisheries  research  vessel,  they  just  happened  to  record  some  fortuitous 
dolphin  signals.  These  signals  were  detected  by  a  sector  side-scanning  sonar  of  high 
bearing  and  time  resolution,  used  as  a  passive  listening  device.  The  acoustic  emissions 
from  the  dolphins  had  significant  energy  at  frequencies  around  305  kHz.  Again,  this  was 
about  one  octave  higher  than  previously  observed.  [Ref.  12] 

C.  RECENT  BROADBAND  MEASUREMENTS  BY  ARL-UT 

The  hydrophone  frequency  response  of  the  prior  measurements  of  Evans,  Au,  and 
Mitson  never  extended  high  enough  to  conclusively  capture  all  of  the  high  frequency 
components.  The  high  resolution  capability  of  cetacean  sonars  prompted  ARL-UT  to 
conduct  further  research  into  the  existence  of  higher  frequencies  that  may  have  been 
overlooked  in  prior  research.  The  "Uncertainty  Principle"  hypothesis  was  proposed  to 
the  Office  of  Naval  Research,  who  then  funded  its  testing  by  the  scientific  method.  The 
two  species  of  dolphins  recorded  in  1997  by  Muir,  Blackwood,  and  Wilson  in  San  Diego 
Bay  were  found  to  emit  significant  high  frequency  signal  components  extending  to  as 
high  as  400  to  500  kHz  [Ref.  1].  These  signals  were  recorded  using  a  hydrophone 
capable  of  measuring  biosonar  signals  up  to  2  MHz.  Details  of  the  hydrophone, 
experimental  configuration,  procedures,  and  results,  are  described  below. 

1.  Wide  Band  Hydrophone  Characteristics 

The  wide  band  hydrophone,  designed  by  Mr.  Lew  Thompson  at  the  Applied 
Research  Laboratories,  University  of  Texas  at  Austin,  was  made  from  a  one  centimeter 
diameter  thin  disk  of  piezo-composite  material.  This  disk  consisted  of  a  mixture  of 
piezo-ceramic  and  a  plastic  material  that  is  inherently  wide  band.  The  transducer 
housing  was  made  of  a  soft,  thin  polyurethane  material.  A  castor  oil  bath  coupling 
medium  was  used  within  the  housing.  Figure  2.4  shows  a  diagram  of  the  hydrophone 
and  the  frequency  response  curve.  Notice  that  the  hydrophone  is  useful  up  to  2  MHz. 
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dB  re  IV/pPa 


There  is  a  notch  present  at  550  kHz,  which  will  be  eliminated  in  future  designs.  Below 
500  kHz,  the  response  curve  is  fairly  flat  and  deviations  from  this  flatness  were  corrected 
for  in  the  data  analysis.  [Ref.  1] 


Figure  2.4  Wide  Band  Hydrophone  and  Receive  Sensitivity  Curve  From 
Ref.[l] 
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2.  Experimental  Configuration 


Recordings  were  conducted  of  the  sounds  emitted  by  two  species  of  captive 
research  dolphins  in  May  and  October  1997,  at  the  Space  and  Naval  Warfare  Systems 
Center  (SPAWARSYSCEN),  in  San  Diego.  The  measurements  utilized  a  bite  bar,  and 
targets  consisting  of  hollow  metal  spheres,  and  bags  of  both  rock  and  junk  metal,  as  seen 
generically  in  Figure  2.5.  First  the  dolphins  were  trained  blindfolded  and  rewarded  fish 
to  eat  for  correctly  identifying  the  different  targets  when  they  were  lowered  in  the  water. 
The  animals  indicated  a  positive  classification  by  emitting  a  whistle,  which  can  be  heard 
by  the  trainers  and  scientists.  Many  data  sets  were  acquired  on  two  bottlenosed  dolphins 
named  Bertha  and  Slooper,  as  well  as  a  beluga  whale  named  Muk  Tuk.  A  random 
sequence  of  designated  real  and  false  targets  were  serially  offered  to  each  animal  during 
the  course  of  the  experiment,  in  order  to  keep  the  animals  alert,  functioning  to  their  best 
capability,  and  to  keep  them  “honest”.  [Ref.  1] 


3.  Experimental  Results 

Experiments  on  the  dolphin  were  conducted  with  the  bite  bar  at  a  depth  of  0.71m, 
with  the  hydrophone  on  axis  with  the  rostrum,  and  located  at  a  range  of  4.65m  from  the 
rostrum.  The  targets  were  placed  at  a  distance  of  9.0m  from  the  bite  bar.  When  the 
target  is  presented  to  the  dolphin,  it  usually  begins  pinging  on  the  target  with  a  series  of 
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rapid  fire  “clicks”,  often  called  a  “click  train.”  A  raw  data  recording  of  a  click  train  is 
shown  in  Figure  2.6.  Notice  that  the  click  train  consists  of  an  increasing  amplitude, 
followed  by  several  pulses  emitted  at  maximum  amplitude,  and  finally  a  decaying 
amplitude.  This  figure  also  shows  a  gradually  increasing  time  between  clicks.  In  this 
example,  there  are  43  clicks  in  a  time  span  of  1.2  seconds.  The  dolphin  clicks  are 
separated  by  about  23  msec,  which,  at  the  speed  of  sound  in  water,  is  the  two  way  travel 
time  to  a  target  located  at  a  range  of  about  21  meters.  For  these  measurements,  the  range 
to  the  target  was  9  meters,  indicating  that  the  echo  from  one  click  was  received  prior  to 
the  emission  of  the  subsequent  click  [Ref.  1]. 


Dolphin  Click  Train 


Figure  2.6  Bottlenose  Dolphin  Echolocation  Click  Train 

Figure  2.7  shows  one  dolphin  click  from  the  middle  of  the  “click  train”  in  Figure 
2.6.  This  raw  data  time  series  contains  128  digitization  points  and  shows  the  relative 
amplitude  of  the  signal  versus  time.  As  can  be  seen,  the  signal  contains  a  few  cycles  of 
powerful  transmission,  followed  by  a  few  weak  cycles  of  acoustic  radiation.  Kamminga 
suggests  that  these  weak  cycles  may  indicate  where  reverberations,  possibly  due  to 
reflections  inside  the  dolphin's  head,  against  the  skull  or  air  sacs,  interfere  with  the  actual 
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first-emitted  sonar  signal  [Ref.  13].  This  time  domain  of  the  dolphin  click  also  shows 
very  sharp  maximum  and  minimum  peaks,  which  may  contain  very  high  frequency 
components. 


Slooper  Dolphin  Click 


Figure  2.7  Bottlenose  Dolphin  Echolocation  Click 


A  proper  acoustic  convention  is  to  show  the  energy  flux  spectral  density  in  a 
signal  with  units  of  pico-joules  per  meter  squared-hertz  [Ref.  14].  The  energy  flux 
spectral  density  of  this  bottlenose  dolphin  click  was  calculated  using  the  MATLAB 
Program  in  Appendix  A,  and  is  presented  in  Figure  2.8.  The  program  accounted  for  the 
frequency  dependence  of  the  ARL-430  hydrophone  receive  sensitivity.  This  figure 
compares  the  click  plus  noise  to  just  noise  (noise  obtained  just  prior  to  the  click  train), 
and  clearly  shows  significant  energy  present  in  the  dolphin  sonar  signal  at  frequencies 
well  above  200  kHz,  much  higher  than  previously  reported.  Utilizing  the  same 
MATLAB  program,  the  energy  flux  spectral  density  was  plotted  for  three  separate  clicks 
from  the  “click  train”  in  Figure  2.6.  These  clicks  were  extracted  from  the  beginning, 
middle,  and  end  of  the  “click  train.”  Figure  2.9  shows  that  the  high  frequency  energy 
flux  increases  near  the  peak  amplitude  of  the  click  train  emissions. 
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Energy  Flux  Spectral  Density  Level  dB  re  [pJ/(nrr  Hz)]  ‘  Energy  Flux  Spectral  Density  Level  dB  re  [pJ/(nrr  Hz)] 


Slooper  Click  Energy  Flux  Spectral  Density  Level  ve  Frequency 


Energy  Flux  Spectral  Density  Level  vs  Frequency 
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Further  analysis  conducted  on  the  dolphin  “click  train”  in  Figure  2.6  shows  a 
great  degree  of  repeatability  among  a  sequence  of  clicks.  Taking  five  consecutive  signals 
from  the  middle  of  the  “click  train”  and  performing  a  cross-correlation,  using  the 
MATLAB  program  in  Appendix  B,  produced  maximum  correlation  coefficients  ranging 
from  96-99%.  The  maximum  correlation  coefficients  of  the  five  signals  can  be  seen  in 
Table  2.1.  Figure  2.10  shows  a  plot  of  the  five  signals  aligned  to  a  position  where  the 
signals  best  correlate  to  each  other.  There  is  little  difference  in  the  time  series  between 
the  five  clicks. 
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Table  2.1  Maximum  Correlation  Coefficients  of  Five  Consecutive  Dolphin  Clicks 


Comparison  of  Five  Consecutive  Dolphin  Clicks 
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This  plot  suggests  that,  if  the  high  frequencies  above  200  kHz  are  present  in  one 
click,  then  the  high  frequencies  must  also  be  present  in  other  clicks  from  the  maximum 
amplitude  portion  of  the  click  train. 

This  chapter  has  briefly  explained  the  dolphin  echolocation  system  and  described 
some  past  research  conducted  on  high  frequency  echolocation  components.  It  has  also 
shown  that  much  higher  frequencies  components  are  in  fact  present  within  a  dolphin 
click.  The  next  chapter  will  explain  the  selection  of  an  absorbing  screen  that  is  being 
used  to  answer  the  question:  is  the  dolphin  using  these  high  frequencies  during 
echolocation? 
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III.  SELECTION  OF  ABSORPTIVE  ACOUSTIC  SCREENS 


The  question  remains  as  to  whether  or  not  the  dolphin  is  actually  using  the  newly 
discovered  high  frequency  components  described  in  Chapter  Two.  In  order  to  test  this 
hypothesis,  Dr.  Muir  proposed  to  ONR  that  an  experiment  be  performed  in  which  an 
acoustic  filter  to  suppress  the  high  frequencies  is  placed  between  a  dolphin  and  a  target, 
and  it  is  determined  whether  the  dolphin's  ability  to  detect  a  target  is  impaired.  This 
chapter  will  describe  the  theory,  selection,  and  laboratory  experiments  conducted  on 
various  materials  for  this  purpose. 

A.  INSERTION  LOSS  AND  REFLECTION  LOSS  THEORY 

Anechoic  coatings,  and  bulk  absorbing  materials  ideally  reflect  zero  percent  of 
sound  incident  upon  them,  and  therefore  present  themselves  as  a  good  absorber.  These 
materials  can  be  evaluated  by  measuring  the  percentage  of  sound  transmitted  through  and 
reflected  from  the  material  when  a  sample  is  immersed  in  water.  These  two  measured 
characteristics,  called  the  “insertion  loss”  and  “reflection  loss,”  are  defined  by  Equations 
3.1,  and  3.2  respectively. 


f 

Insertion  loss  =  20  log 

V 


Incident  rms  sound  pressure  N 
Transmitted  rms  sound  pressure , 


Reflection  loss  =  20  log 


Incident  rms  sound  pressure 
Reflected  rms  sound  pressure 


(3.1) 


(3.2) 


Plane-wave  propagation  is  assumed  in  both  definitions,  and  both  characteristics  are 
expressed  in  positive  decibel  units.  [Ref.  15] 

Insertion  loss  is  the  reduction  in  the  signal,  in  decibels,  caused  by  inserting  the 
material  between  the  sound  source  and  the  receiver,  with  diffraction  and  refraction  effects 
absent.  The  insertion  loss  of  the  material  is  due  to  the  combination  of  sound  reflected 
from  the  material  and  sound  absorbed  in  the  material,  as  illustrated  in  Figure  3.1 .  It  is 
also  a  function  of  the  material's  dimensions  and  properties.  When  plane  waves  impinge 
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Figure  3.1  Reflection  and  Transmission  of  Plane  Waves  Normally  Incident  on  a 
Layer  After  Ref.  [3]. 


normally  on  a  plate  of  a  homogeneous  absorbing  material  of  uniform  thickness,  and  water 
is  on  both  sides  of  the  plate,  the  theoretical  insertion  loss  is  given  by: 


IL  =  20 log  (w  + 1  ~  Jlf  exPO'2fc^) -(m-\  +  jrf  exp(-2ozQ  (3  3) 

4m(l  -  jr)  exp (jkd  -  ad) 


where  m  is  the  real  ratio  of  the  characteristic  impedance  of  the  material  to  the 
characteristic  impedance  of  water  ( pc/p0c0 ),  k  is  the  wave  number  (co/c),  d  is  the  thickness 
of  the  sample  material  in  meters,  a  is  the  longitudinal  attenuation  constant  in 
Nepers/meter,  and  r  is  the  loss  parameter  of  the  material,  defined  by  Equation  3.4. 


In  applying  Equation  3.3  in  this  present  investigation  to  model  the  behavior  of  the  sound 
absorbing  materials,  a  was  taken  to  be  at  most  a  quadratic  function  of  the  frequency  f: 

a{f)  =  af  +  bf 2  (3.5) 
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The  coefficients  a  and  b  can  be  determined  by  curve  fitting  Equation  3.3  to  the 
experimental  data  for  the  acoustic  absorbing  panel,  or  they  can  be  measured  directly  in  a 
large  material  sample.  [Ref.  1 6] 

Reflection  loss  is  a  measure  of  how  much  the  incident  sound  pressure  level  is 
reduced  upon  reflection,  and  depends  on  the  acoustical  impedance  mismatch  at  the 
reflection  boundary.  The  boundary  acoustical  impedance  in  turn  depends  on  the  material 
itself,  and  its  thickness.  Anechoic  coatings  often  consist  of  a  layer  of  rubber  material 
mixed  with  either  low-acoustic-impedance  air  voids  or  high-acoustic-impedance  metal 
particles  or  both.  The  coating  should  have  a  good  impedance  match  with  the  water  so  the 
sound  energy  will  not  be  reflected,  and  after  entering  the  coating,  the  sound  energy 
should  optimally  be  absorbed,  otherwise  the  sound  will  be  reflected  at  some  subsequent 
boundary.  [Ref.  15] 

Theoretical  reflection  loss  can  be  computed  for  sound  incident  normally  on  a 
uniform  plate  or  panel  of  absorbent  material  immersed  in  water,  and  is  given  by: 


R£=201og 


[(ot  +  1-  jr)l(m-\+ jr)]exp(Jkd+  cd)-[(m-\+ jr)/(m  +  \-jr)]exp(-jkd-ad) 
exp  (jkd+  cd)  -  expf- jkd- ad) 


(3.6) 


where  m,  kd,  a,  and  r  are  the  same  as  in  Equation  3.3.  In  applying  Equation  3.6  to  model 
the  behavior  of  the  sound  absorbing  material,  a  was  again  taken  to  be  a  quadratic 
function  of  frequency.  Reflection  loss  is  a  maximum  for  frequencies  at  which  the  panel 
thickness  is  a  multiple  of  a  half  wavelength.  [Ref.  16] 


B.  LIQUID  MOLDING  COMPOUND  INVESTIGATION 


Attempting  to  find  a  material  with  the  correct  properties  to  act  as  a  low  pass  filter 
proved  quite  challenging.  Experimentation  first  involved  the  testing  of  over  the  counter 
liquid  molding  compounds.  The  liquid  molding  was  blended  in  a  household  blender  to 
introduce  air  bubbles.  The  molding  compound  was  then  poured  into  a  wooden  mold  and 
allowed  to  set  for  12  hours.  As  the  material  solidified,  the  air  bubbles  migrated  toward 
the  top  of  the  material.  The  final  result  was  a  substance  composed  of  small  air  bubbles  at 
the  bottom  of  the  sample  and  larger  bubbles  near  the  top.  This  material  was  expected  to 
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act  as  a  lossy  acoustic  medium  allowing  sound  to  transmit  through  the  small  bubble  side 
with  attenuation,  and  reflecting  off  the  large  bubble  side.  Measuring  the  insertion  loss 
on  numerous  tiles  with  different  thickness  and  bubble  consistency  showed  the  material  to 
be  nonabsorptive  and  purely  reflective  over  the  frequency  range  of  60  to  400  kHz. 
Insertion  loss  experiments  were  also  conducted  on  open  cell  foam  as  recommend  by 
Undersea  Warfare  Center,  Newport  RI.  Again  these  materials  were  purely  reflective  and 
non-absorbing  over  the  same  frequency  range.  The  search  for  an  ideal  absorbing  tile 
ended  when  two  pieces  of  SOAB,  sound  absorbing  material ,  were  discovered  in  one  of 
the  acoustic  laboratories  at  the  Naval  Postgraduate  School. 

C.  SOAB  PROPERTIES  AND  TESTING 

Sound  absorbing  material  (SOAB)  was  first  invented  by  the  Germans  in  World 
War  II.  The  German  "Alberich"  coating,  consisting  of  a  rubber  layer  with  air-filled 
voids,  was  experimentally  cemented  to  the  outside  of  a  number  of  U-boats  [Ref.  17]. 
SOAB  is  a  porous  panel  made  from  butyl  rubber  imbedded  with  aluminum  powder. 
Microscopic  air  bubbles  become  attached  to  the  aluminum  powder  surface  in  the 
manufacturing  process.  The  sound  absorption  is  dependent  upon  these  air  bubbles.  In 
the  1950's,  B.  F.  Goodrich  Company  of  Akron,  Ohio,  made  commercially  available  three 
different  types  of  sound  absorbing  panels:  SOAB  I,  SOAB  II,  and  SOAB  III.  Each  of 
these  coatings  consisted  of  different  aluminum  powder  loads  which  varied  the  material's 
density.  The  absorbing  material  was  researched  as  anechoic  coatings  for  acoustic  test 
tank  lining  to  absorb  reflections.  In  1961,  B.  F.  Goodrich  produced  a  plot  showing  the 
transmission  loss  versus  frequency  for  a  SOAB  baffle  as  a  function  of  thickness.  The 
plot,  shown  in  Figure  3.2,  did  not  mention  the  type  of  SOAB,  and  the  frequency  range 
only  extended  up  to  60  kHz.  [Ref.  1 8] 

The  first  important  parameter  measured  in  order  to  determine  the  effectiveness  of 
SOAB  for  use  as  an  acoustic  filter  was  the  insertion  loss.  Tests  were  performed  in  a 
small  anechoic  tank  designed  for  ultrasonic  research.  The  four  sides  and  bottom  of  the 
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Figure  3.2  SOAB  Transmission  Loss  versus  Frequency  From  Ref.  [18]. 

tank  are  shaped  like  “butterfly  wings”  to  reduce  the  effect  of  standing  waves  and 
reflections,  as  described  by  Willette  and  Muir  [Ref.  19].  The  experiment  consisted  of  a 
transducer  projecting  a  tone  burst  of  sound  toward  a  hydrophone,  and  measuring  the 
received  hydrophone  output  voltage.  The  experiment  geometry  can  be  seen  in  Figure  3.3. 
The  hydrophone  was  placed  30  cm  from  the  transducer  to  ensure  meeting  far  field 
criteria.  The  absorbing  filter  was  placed  between  the  projector  and  the  receiver  and  the 
voltage  was  again  recorded.  The  insertion  loss  was  calculated  by  substituting  the 
measured  voltages  for  the  pressure  terms  and  rewriting  Equation  3.1  as: 


21 


(3.7) 


7Z,  =  201og— 

V, 

where  IL  is  the  insertion  loss  in  decibels,  V,-  is  the  direct  hydrophone  receive  voltage,  and 
Vt  is  the  hydrophone  receive  voltage  through  the  material. 


Figure  3.3  Geometry  For  Insertion  Loss  Experiment 


The  second  parameter  measured  for  the  acoustic  filter  was  reflection  loss. 
Minimizing  reflections  off  the  absorbing  filter  is  desirable,  so  as  not  to  contuse  the 
dolphin  with  reflected  signals  other  than  from  the  target.  These  measurements  were 
performed  in  the  same  anechoic  tank  as  described  above.  The  experiment  consisted  of  a 
transducer  projecting  pulsed  sound,  a  hydrophone,  and  a  sample  tile  as  shown  in  Figure 
3.4.  The  hydrophone  was  placed  50cm  from  the  projector,  without  the  filter  present,  and 
the  receive  voltage  was  recorded.  The  filter  was  then  placed  in  the  tank  at  50  cm  from 
the  projector  and  the  hydrophone  moved  to  25cm  ,  halfway  between  the  projector  and  the 
filter.  The  hydrophone  receive  voltage  was  again  recorded.  Correcting  for  spherical 
spreading,  the  reflection  loss  was  calculated  using: 

RL  =  201og  p-  (3.8) 

Y  r 

where  RL  is  the  echo  reduction  in  decibels,  and  V,  and  Vr  are  the  hydrophone  receive 
voltages  with  and  without  the  filter  respectively. 
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Figure  3.4  Geometry  For  Reflection  Loss  Experiment 


D.  EXPERIMENTAL  RESULTS 

The  results  from  the  insertion  loss  experiment  are  shown  in  Figure  3.5.  Insertion 
Loss  was  measured  on  SOAB  I  ‘//'-thick  sample,  SOAB  III  ‘/-"-thick  sample,  and  for 
both  SOAB  I  and  SOAB  III  together,  making  a  %"-thick  sample.  Figure  3.5  shows  that 
increasing  the  thickness  of  the  SOAB  screen  also  increases  the  insertion  loss.  As  stated 
in  the  beginning  of  this  chapter,  it  is  desired  to  use  these  absorbing  filters  to  suppress  the 
high  frequencies  above  100-150  kHz  and  determine  the  dolphin's  ability  to  detect  targets 
with  and  without  the  filters.  Of  the  choices  available,  the  optimal  filter  to  use  in  this 
application  would  be  the  SOAB  III  ‘/"-thick,  because  it  has  the  least  amount  of 
absorption  at  the  lower  frequencies,  and  yet  provides  a  greater  absorption  at  the  higher 
frequencies  to  suppress  them.  It  should  be  mentioned  that  the  total  insertion  loss  is  twice 
that  shown  in  Figure  3.5  due  to  the  two  way  travel  through  the  absorbing  screen. 

Since  the  absorbing  screens  would  be  hanging  underwater  during  the  actual 
experiment  with  the  dolphin,  it  was  of  interest  to  measure  the  insertion  loss  as  a  function 
of  incidence  angle.  Even  though  the  screen  would  be  tethered  so  as  to  restrict  its 
movement,  there  might  be  slight  movement  during  the  open  ocean  experiment.  Figure 
3.6  shows  the  Insertion  Loss  of  the  ‘/’’-thick  screen  as  a  function  of  incidence  angle  off 
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axis.  As  one  can  see,  there  is  only  a  slight  variation  in  insertion  loss  with  incidence  angle 
for  nearly  normal  incidence. 


50  100  150  200  250  300  350  400 

frequency  (kHz) 

Figure  3.5  Insertion  Loss  vs  Frequency  For  Three  SOAB  Samples 


SOAB  Insertion  Loss  As  A  Function  Of  Incidence  Angle  Off  Axis 


Figure  3.6  Insertion  Loss  vs.  Incidence  Angle  Off  Axis  For  V”  SOAB 
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In  comparing  theory  (Equations.  3.3  and  3.6)  to  recorded  measurements,  it  was 
necessary  to  obtain  an  accurate  value  for  the  ratio  of  acoustic  impendances  m,  and  to 
measure  the  impedance  for  both  pieces  of  absorbing  filters.  First,  the  density  of  each 
material  was  calculated  by  weighing  a  small  sample  of  the  material  and  then  submersing 
the  piece  in  water,  and  measuring  the  volume  of  water  displaced.  The  weight  divided  by 
volume  of  SOAB  I  (14")  and  SOAB  III  (14")  resulted  in  densities  of  1520  kg/m3,  and  993 
kg/m3,  respectively. 

The  speed  of  sound  in  each  material  was  estimated  from  measurements  of  the 
time  change  or  phase  shift  of  the  received  signal  zero  axis  crossings  with  and  without  the 
filter  present.  The  filter  was  placed  directly  in  front  of  the  hydrophone,  with  its  shortest 
axis  being  parallel  to  the  direction  of  propagation.  The  time  change  of  the  zero  axis 
crossings  was  measured  near  the  beginning  and  middle  of  the  received  tone  burst  at  60 
kHz,  and  130  kHz,  with  an  accuracy  of  plus  or  minus  0.01  microseconds.  The  filter 
sound  velocity  could  then  be  calculated  by: 

C  C-A  (3.9) 

'  ds  -  CwAt 

where  Cs  and  Cw  are  the  sound  speed  in  the  filter  and  water  respectively,  ds  is  the  filters 
thickness,  and  At  is  the  time  shift.  Using  Equation  3.9  and  the  averaged  time  shift 
produced  a  velocity  of  1586  m/s  for  SOAB  I,  and  1662  m/s  for  SOAB  III. 

The  measured  insertion  loss  for  the  14”  SOAB  filter  is  plotted  against  theory  in 
Figure  3.7.  The  theory  curve  was  calculated  using  Equation  3.3,  and  making  a  number  of 
successive  approximations  for  the  coefficients  a  and  b  in  Equation  3.5,  until  the  theory 
best  matched  the  experimental  data  visually.  Theory  and  experimental  matched  best 
when  values  of  a  and  b  where  chosen  as  0.20  Np/m/kHz2,  and  12  x  10‘5  Np/m/kHz2, 
respectively. 

It  was  also  necessary  to  measure  reflection  loss  in  our  14”  SOAB  filter,  to 
minimize  the  reflected  dolphin  echolocation  signal  from  the  front  surface  of  the  filter. 
Reflection  Loss  was  measured  and  compared  with  theory  in  Figure  3.8.  The  theory  curve 
was  calculated  using  Equation  3.6,  and  again  choosing  the  same  values  for  a  and  b  in 
equation  3.5,  until  theory  best  matched  experimental  data  visually.  This  figure  clearly 
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shows  that  reflection  loss  theory  is  a  maximum  for  frequencies  at  which  the  absorptive 
screen  thickness  is  a  multiple  of  a  half  wavelength.  The  reflection  loss  shows  that  the 
screen  is  non-reflecting,  even  at  the  higher  frequencies. 


Insertion  Loss  (1/2"  SOAB  Hie) 


Figure  3.7  Experimental  And  Theory  Insertion  Loss 


Reflection  Loss  (1/2"  SOAB  Tile) 


Figure  3.8  Experimental  and  Theory  Reflection  Loss 
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This  chapter  has  explained  insertion  loss  and  reflection  loss  theory  and  applied  it 
to  experimental  data.  It  has  also  shown  that  sound  absorbing  tile  (SOAB)  is  the  best 
material  for  use  as  a  low  pass  filter.  The  next  chapter  will  explain  the  procedure  in  using 
this  screen  to  filter  out  the  dolphins  echolocation  high  frequencies  and  the  effect  of  its 
insertion  on  the  animals  target  detection  performance. 
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IV.  BOTTLENOSE  DOLPHIN  ECHOLOCATION  EXPERIMENT 


This  chapter  will  discuss  the  configuration,  procedure,  and  results  from  the 
bottlenose  dolphin  echolocation  experiment  performed  in  July  1998  at  the  Space  and 
Naval  Warfare  Systems  Center  (SPAWARSYSCEN),  in  San  Diego.  The  experiment 
involved  inserting  pieces  of  sound  absorbing  filters  (SOAB)  in  front  of  the  dolphin  and 
observing  the  effect  on  the  animal's  target  detection  performance. 

A.  EXPERIMENTAL  CONFIGURATION 

The  experiment,  as  shown  in  Figure  4.1,  was  performed  in  a  floating  pen, 
surrounded  by  fish  netting  and  open  to  San  Diego  Harbor.  The  bite  plate  support,  ARL- 
430  hydrophone,  and  the  SOAB  filters  were  attached  on  one  side  of  the  pen.  The  targets 
were  hung  from  the  other  side  of  the  pen.  The  bite  plate  support,  as  shown  in  Figure  4.2, 
was  constructed  from  2”  PVC  pipes  and  contained  a  sliding  neoprene  door.  When  the 
neoprene  door  was  lowered  by  the  trainer,  there  was  provided  an  unobstructed  aperture 
through  which  the  dolphin  could  project  its  sonar  signals.  When  raised,  the  neoprene 
door  blocked  the  animal’s  sonar  signals  from  the  hydrophone  and  targets.  The  bite  plate 
was  fixed  at  a  depth  of  78  cm. 


29 


Figure  4.2  Photo  of  Bite  Plate  And  Neoprene  Door  Assembly 

The  ARL-430  hydrophone,  shown  in  Figure  4.3,  was  placed  1.4  meters  from  the 
bite  plate  and  on  axis  at  a  depth  of  78cm.  The  hydrophone  was  connected  to  a  12-bit 
National  Instruments  data  acquisition  card,  which  converted  analog  signals  to  digital 
signals.  The  card  was  then  connected  to  a  personal  computer  for  data  collection. 


Figure  4.3  Photo  of  ARL-430  Hydrophone  And  Preamplifier  Next  to  a  Penny 
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The  SOAB  filters  were  attached  to  square  frames  made  from  1”  PVC  piping.  All 
of  the  PVC  piping  had  holes  drilled  in  them  which  allowed  the  pieces  to  be  filled  with 
seawater  as  they  were  submerged.  The  filters  were  hung  1.8  m  from  the  bite  plate  and 
the  middle  of  the  filter  was  on  axis  at  a  depth  of  78  cm.  The  filters  were  rigged  such  that 
the  down  position  was  well  below  the  bite  plate  support  and  out  of  the  dolphins  sonar 
beam.  In  the  up  position  the  SOAB  filters  were  centered  on  axis  with  the  bite  plate  so 
that  the  animal's  sonar  beam  had  to  penetrate  the  filter. 

The  “target,”  as  shown  in  Figure  4.4,  consisted  of  a  8cm  diameter,  stainless  steel, 
water-filled  sphere.  A  thin  monofilament  line  was  used  to  lower  the  target  to  a  depth  of 
78  cm  or  raise  it  out  of  the  water.  The  “target  absent”  line  contained  just  a  lead  weight 
with  no  steel  sphere.  Lead  weights  were  hung  at  the  same  depth  from  both  lines  to  keep 
tension  on  the  lines.  The  lead  weight  on  the  “target  present”  line  was  attached  a  meter 
below  the  steel  ball.  Both  targets  were  hung  at  a  distance  of  7  m  from  the  bite  plate.  A 
photo  of  the  experimental  configuration  is  shown  in  Figures  4.5. 


Figure  4.4  Photo  of  Stainless  Steel  Water  Filled  Sphere  Target 
Next  To  A  4x4  Wood  Post 
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Figure  4.5  Photo  of  Experiment  Configuration  Showing  Bite  Plate  and  Target  Stand 


B.  EXPERIMENTAL  PROCEDURE 

Figure  4.6  shows  the  trainer  first  placing  suction  cups  over  the  dolphin’s  eyes  to 
prevent  it  from  seeing  during  the  trials.  The  dolphin  was  then  required  to  station  on  the 
bite  plate  while  the  neoprene  door  remained  in  the  up  position.  A  trial  started  when  the 
trainer  lowered  the  door,  cueing  the  dolphin  to  commence  its  sonar  search.  After 
completing  the  sonar  search,  the  dolphin  would  respond  with  a  whistle  to  indicate  “target 
present”  response,  or  remain  quiet  to  indicate  “target  absent”  response.  If  the  dolphin 
provided  the  correct  response,  the  trainer  would  signal  the  dolphin  with  a  “bridge.”  A 
"bridge"  is  a  high  pitch  whistle  that  signals  the  dolphin  to  leave  the  bite  plate  and  return 
to  the  surface  for  a  reward  (usually  fish  or  squid).  If  the  dolphin  gave  the  incorrect 
response,  the  trainer  would  signal  a  “delta”  to  the  dolphin.  A  "delta"  is  a  tone  that 
informs  the  dolphin  it  made  the  wrong  choice  and  would  not  receive  any  reward.  Figure 
4.7  shows  the  trainer  signaling  the  dolphin  with  a  device,  that  produces  the  high  pitch 
whistle  or  tone,  known  as  a  “Bundy  Box.” 

Four  sets  of  20  trials  were  performed  on  one  bottlenose  dolphin  named  Slooper. 

A  random  sequence  of  trials  presenting  combinations  of  target,  or  no  target,  filter,  or  no 


filter,  was  serially  offered  to  the  animal.  The  first  set  of  20  trials  consisted  of  an  empty 
frame  with  no  SOAB  filter.  The  second  set  presented  a  framed  54”  SOAB  filter.  The 
third  set  presented  a  framed  54”  SOAB  filter.  The  last  set  presented  a  framed  %”  SOAB 
filter  which  was  constructed  from  the  54”  and  54”  SOAB  tiles  mounted  on  one  frame. 


Figure  4.7  Photo  of  Trainer  Signaling  Dolphin  With  Bundy  Box 
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c. 


EXPERIMENTAL  RESULTS 


1.  Target  Recognition 


The  results  from  the  three  sets  of  trials  containing  the  framed  SOAB  filters  are 
shown  in  Tables  4. 1  through  4.3.  Incorrect  responses  are  shown  highlighted.  Looking 
at  the  results  from  the  three  sets  of  data,  one  might  first  conclude  that  the  *4”  SOAB  filter 
produced  little  effect  to  the  dolphin’s  echolocation  ability,  while  the  14”  and  %”  SOAB 
filters  possibly  impaired  the  dolphin’s  echolocation  ability.  A  closer  look  at  the  data 
reveals  that  the  SOAB  filters,  in  this  experimental  configuration,  probably  had  little  or  no 
effect  on  the  dolphin's  echolocation  ability.  In  the  first  test,  using  the  14”  SOAB  filter, 
the  dolphin  was  perfect  in  identifying  the  presence  of  a  target  from  no  target.  In  the 
second  test,  using  the  14”  SOAB  filter,  the  dolphin  was  incorrect  in  four  responses.  Three 
of  these  four  incorrect  responses  occurred  with  no  filter  present,  and  possibly  was  a  result 
of  the  animal  being  distracted  during  the  individual  trial.  In  the  third  test,  using  the  %” 
SOAB  filter,  the  dolphin  gave  only  one  incorrect  response,  and  again  occurred  when  no 
filter  was  present. 


TRIAL 

TARGET 

FILTER 

DOLPHIN 

DOLPHIN 

NUMBER 

PRESENT 

PRESENT 

RESPONSE 

DECISION 

1 

YES 

NO 

WHISTLE 

CORRECT 

2 

YES 

NO 

WHISTLE 

CORRECT 

3 

YES 

NO 

WHISTLE 

CORRECT 

4 

NO 

YES 

QUIET 

CORRECT 

5 

NO 

YES 

QUIET 

CORRECT 

6 

NO 

YES 

QUIET 

CORRECT 

7 

NO 

NO 

QUIET 

CORRECT 

8 

YES 

YES 

WHISTLE 

CORRECT 

9 

NO 

NO 

QUIET 

CORRECT 

10 

NO 

NO 

QUIET 

CORRECT 

11 

NO 

YES 

QUIET 

CORRECT 

12 

YES 

YES 

WHISTLE 

CORRECT 

13 

NO 

YES 

QUIET 

CORRECT 

14 

YES 

NO 

WHISTLE 

CORRECT 

15 

YES 

YES 

WHISTLE 

CORRECT 

16 

YES 

YES 

WHISTLE 

CORRECT 

17 

YES 

NO 

WHISTLE 

CORRECT 

18 

NO 

NO 

QUIET 

CORRECT 

19 

NO 

NO 

QUIET 

CORRECT 

20 

YES 

YES 

WHISTLE 

CORRECT 

20  -  Correct  Responses  0  —  Incorrect  Responses 

Table  4.1  Echolocation  Test  Using  14"  SOAB  Tile 
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TRIAL 

TARGET 

FILTER 

DOLPHIN 

DOLPHIN 

NUMBER 

PRESENT 

PRESENT 

RESPONSE 

DECISION 

1 

YES 

NO 

WHISTLE 

CORRECT 

2 

YES 

NO 

OUIET 

WRONG 

3 

YES 

NO 

WHISTLE 

CORRECT 

4 

YES 

YES 

OUIET 

WRONG 

5 

NO 

YES 

OUIET 

CORRECT 

6 

NO 

YES 

OUIET 

CORRECT 

7 

NO 

NO 

WHISTLE 

WRONG 

s 

NO 

YES 

OUIET 

CORRECT 

9 

NO 

NO 

OUIET 

CORRECT 

10 

YES 

YES 

WHISTLE 

CORRECT 

11 

NO 

YES 

OUIET 

CORRECT 

12 

YES 

YES 

WHISTLE 

CORRECT 

13 

YES 

NO 

WHISTLE 

CORRECT 

14 

NO 

NO 

WHISTLE 

WRONG 

15 

YES 

YES 

WHISTLE 

CORRECT 

16 

NO 

NO 

OUIET 

CORRECT 

17 

YES 

NO 

WHISTLE 

CORRECT 

18 

YES 

YES 

WHISTLE 

CORRECT 

19 

NO 

NO 

OUIET 

CORRECT 

20 

NO 

YES 

OUIET 

CORRECT 

16-  Correct  Responses  4  -  Incorrect  Responses:  2  false  positives,  and 

2  false  negatives 

Table  4.2  Echolocation  Test  Using  V”  Filter 


FILTER 

PRESENT 

DOLPHIN 

RESPONSE 

DOLPHIN 

DECISION 

i 

YES 

NO 

WHISTLE 

CORRECT 

2 

NO 

NO 

OUIET 

CORRECT 

3 

YES 

YES 

WHISTLE 

CORRECT 

4 

NO 

NO 

OUIET 

CORRECT 

5 

YES 

NO 

WHISTLE 

CORRECT 

6 

NO 

YES 

OUIET 

CORRECT 

7 

YES 

NO 

WHISTLE 

CORRECT 

8 

NO 

NO 

OUIET 

CORRECT 

9 

YES 

YES 

WHISTLE 

CORRECT 

10 

YES 

NO 

WHISTLE 

CORRECT 

11 

NO 

YES 

OUIET 

CORRECT 

12 

NO 

YES 

OUIET 

CORRECT 

13 

YES 

NO 

OUIET 

WRONG 

14 

YES 

YES 

WHISTLE 

CORRECT 

15 

NO 

NO 

OUIET 

CORRECT 

16 

NO 

NO 

OUIET 

CORRECT 

17 

NO 

YES 

OUIET 

CORRECT 

18 

YES 

YES 

WHISTLE 

CORRECT 

19 

NO 

YES 

OUIET 

CORRECT 

20 

YES 

YES 

WHISTLE 

CORRECT 

19  -  Correct  Responses  1  -  Incorrect  Response:  false  negative 

Table  4.3  Echolocation  Test  Using  V”  Filter 
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2. 


Time  Duration  Of  Click  Trains 


The  duration  of  the  dolphin  click  train  (average  time),  in  seconds,  is  shown  in 
Table  4.4.  The  rows  in  the  table  represent  each  of  the  recorded  data  sets,  where  a 
different  size  filter  was  attached  to  the  PVC  frame.  The  columns  represent  the  different 
presentations  of  target  and  filter  present  or  not  present  during  the  individual  data  set  of  20 
trials.  Comparing  the  first  two  columns,  the  click  train  is  shorter  when  the  target  and 
filter  were  both  presented  to  the  dolphin.  This  indicates  that  the  dolphin  was  able  to 
identify  the  presence  of  the  target  just  the  same  with  or  without  the  filter  present.  When 
there  was  no  target  present,  the  click  train  length  was  considerably  longer  indicating  that 
the  animal  was  really  searching  for  the  target  before  it  decided  it  was  absent.  The  last 
two  columns  show  that  the  click  train  length  increased  even  more  when  the  filter  and  no 
target  was  presented  to  the  dolphin.  This  trend  appears  to  suggest  that  the  filter  caused 
some  difficulty  in  the  dolphin's  ability  to  ensure  that  the  target  was  in  fact  absent. 


DATA  SET 

TARGET 

NO  FILTER 

TARGET 

FILTER 

NO  TARGET 
NO  FILTER 

NO  TARGET 
FILTER 

EMPTY  FRAME 

1.37  sec 

0.792  sec 

1.29  sec 

1.72  sec 

y*n  FILTER 

1.11  sec 

0.742  sec 

1.76  sec 

2.06  sec 

J4"  FILTER 

1.07  sec 

1 .02  sec 

1.62  sec 

2.18  sec 

%"  FILTER 

1.15  sec 

1.12  sec 

2.14  sec 

2.19  sec 

Table  4.4  Dolphin  Click  Train  Average  Time  During  Echolocation  Trials 


3.  Click  Energy  Flux  Spectral  Density 


Figures  4.8  through  4.1 1  show  eight  clicks  overlapping  from  the  maximum 
amplitude  portion  of  a  click  train  and  its  corresponding  energy  flux  spectral  density.  The 
standard  deviation  of  the  eight  signal's  energy  flux  spectral  density  is  also  shown  on  the 
plot  as  a  dot-dashed  line.  The  signals  analyzed  were  chosen  from  a  trial,  representing 
each  of  the  four  data  sets,  when  a  target  was  present  and  the  dolphin  correctly  identified 
the  target.  The  figures  show  an  increase  in  the  time  series  click  amplitude  and  respective 
energy  flux  spectral  density  as  the  thickness  of  filter  increases.  The  energy  present  at  the 
peak  frequency  around  130  kHz  remains  constant  through  the  four  trials,  but  the  energy 
between  250  kHz  and  500  kHz  increases  as  the  dolphin  attempts  to  penetrate  the  thicker 
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SOAB  filter  with  its  sonar.  It  appears  that  the  dolphin  is  utilizing  a  portion  of  die 
higher  frequency  components  to  identify  the  presence  of  a  target.  The  dolphin  is  also 
having  to  increase  its  signal  amplitude  to  penetrate  the  thicker  filter. 


1 1 _ i - i - 1 - 1 - 1 - 1 

0  20  40  60  80  100  120 
Time  in  Microseconds 


! _ i - i - i - 1— - 1 - 1 

0  100  200  300  400  500  600 


Frequency(kHz) 

Figure  4.8  Eight  Dolphin  Clicks  And  Corresponding  Average  Energy  Flux 
Spectral  Density  With  No  Filter  Present  (±  1  a  shown  as  dot-dash  line) 
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Frequency  (kHz) 


Figure  4.9  Eight  Dolphin  Clicks  And  Corresponding  Average  Energy  Flux 
Spectral  Density  With  14"  Filter  Present  (±  1  cr  shown  as  dot-dash  line) 


Relative  Amplitude 


4.  Diffraction  As  A  Source  Of  Experimental  Uncertainty 

One  possible  explanation  of  why  the  dolphin's  echolocation  performance 
appeared  to  be  unaffected  during  the  insertion  of  acoustic  high-frequency  filters  could  be 
that  a  portion  of  the  dolphin's  sonar  beam  passed  around  the  filter,  as  depicted  in  Figure 
4.12. 


Au  calculated  the  average  3-dB  and  1 0-dB  vertical  beam  widths  of  three 

o  o 

bottlenose  dolphin  transmitted  sonar  beams,  centered  at  130  kHz,  to  be  10.2  and  22.5 
respectively  [Ref.  2].  The  width  of  a  sonar  beam  at  the  distance  of  the  filter  (1.7m)  can 
be  calculated  from  the  relation  s=rO,  where  s  is  the  surface  area,  r  is  the  radius,  and  6  is 
the  beam  width  in  radians.  Using  this  relation,  the  vertical  extent  of  the  dolphin's  sonar 
beam  at  the  filter  would  have  been  30  cm  for  the  3-dB  beam  width,  and  68  cm  for  the 
10-db  beam  width.  The  width  of  the  14"  SOAB  filter  was  53cm,  and  the  width  of  the  14" 
SOAB  filter  was  28cm.  These  results  indicate  that  a  portion  of  the  dolphin's  sonar  beam 
was  not  entirely  blocked  by  the  filter.  No  measurements  were  recorded  in  this 
experiment  to  indicate  whether  or  not  diffracted  paths  were  important.  Future 
experiments  must  ensure  the  filter  is  placed  closer  to  the  dolphin  to  ensure  the  filter 
blocks  the  entire  sonar  beam. 
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This  chapter  has  discussed  the  configuration,  procedure  and  some  significant 
results  from  an  echolocation  target  detection  task  with  the  presence  of  sound  absorbing 
filters.  Although  this  experiment  has  not  provided  conclusive  evidence  on  the  dolphin's 
use  of  high  frequency  components  during  echolocation,  it  has  shown  clear  evidence  that 
some  compensation  was  made  by  the  dolphin  for  the  sound  absorption  caused  by  the 
insertion  of  the  filter. 
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V. 


CONCLUSIONS  AND  RECOMMENDATIONS 


This  final  chapter  briefly  synopsizes  the  noteworthy  observations  made  during  the 
conduct  of  the  research  described  in  this  thesis  and  the  conclusions  that  can  be  drawn 
from  these  observations.  It  also  provides  further  recommendations  for  continuing 
research  efforts. 

The  first  objective  of  this  thesis  was  to  verify  the  existence  of  high  frequency 
components  in  bottlenose  dolphin  echolocation  signals.  The  second  chapter  clearly 
showed  the  existence  of  high  frequency  energy  components  in  a  dolphin  click.  Using  a 
wide-band  hydrophone  allowed  the  high  frequency  components  to  become  visible.  The 
comparison  of  a  sequence  of  clicks  from  the  steady  maximum  amplitude  portion  of  the 
click  train  also  showed  that  if  the  high  frequency  components  are  present  in  one  click, 
then  they  are  present  in  adjacent  clicks,  to  an  amazing  degree  of  uniformity  and 
repeatability. 

The  second  objective  was  to  obtain  and  test  an  absorbing  material  capable  of 
filtering  out  the  higher  frequency  components  in  the  dolphin’s  sonar  signal.  The  third 
chapter  concluded  that  sound  absorbing  material  (SOAB)  best  resembles  a  low  pass  filter 
capable  of  absorbing  a  major  portion  of  the  high  frequency  sonar  signal  components. 

The  material  was  also  found  to  have  very  little  reflectivity.  Acoustically,  this  material 
would  be  ideal  since,  in  theory,  the  dolphin  would  unlikely  be  able  to  detect  the  presence 
of  the  filter  within  its  sonar  beam. 

The  final  objective  was  to  observe  whether  the  dolphin’s  echolocation  ability  was 
impaired  by  the  presence  of  sound  absorbing  material  placed  in  its  beam,  and  so  to 
determine  whether  the  high  frequency  components  are  being  used  by  the  animal  for  sonar 
data  acquisition.  The  fourth  chapter  discussed  the  first  time  a  dolphin  echolocation  target 
detection  experiment  was  conducted  using  various  dimensions  of  SOAB  material  within 
the  dolphin’s  sonar  beam.  Even  though  the  results  of  this  one  experiment  where 
inconclusive  in  answering  the  third  thesis  objective,  it  did  accomplish  filter  testing  with  a 
live  animal  where  follow-on  efforts  can  begin.  It  was  difficult  to  conclude  if  the  dolphin 
was  impaired  at  all  by  the  presence  of  the  various  SOAB  filters.  The  dolphin  was 
confident  in  reporting  the  presence  of  a  target  using  all  three  filters.  The  animal  was  also 
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definitely  compensating  for  the  filters,  by  increasing  its  sound  energy  output,  especially 
for  frequencies  above  100  kHz.  The  dolphin’s  few  errors  may  have  been  a  result  of 
animal  distraction  or  from  being  exposed  to  a  new  experimental  configuration  that  the 
animal  was  not  accustomed  to.  It  is  possible  that  a  small  portion  of  the  animal’s  sonar 
beam  was  not  blocked  by  the  filters,  due  to  their  position  in  the  experiment.  The  only 
way  to  determine  the  true  effect  of  the  filters  on  the  dolphin’s  echolocation  ability  is  to 
conduct  many  more  trials  with  different  experimental  configurations. 

Many  questions  remain  unanswered  in  reference  to  the  dolphin  sonar  system. 
Analysis  of  the  echolocation  data  collected  in  1997  and  recent  1998  trials  could  provide 
even  further  direction  toward  future  research. 

This  thesis  has  provided  a  stepping  stone  toward  future  work  in  discovering  an 
alternative  mine-hunting  system,  which  is  less  expensive  and  has  a  higher  search  rate 
than  the  marine  mammal. 
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APPENDIX  A.  ENERGY  FLUX  SPECTRAL  DENSITY  ANALYSIS 

This  appendix  contains  the  MATLAB  programs  used  to  conduct  the  energy  flux 
spectral  density  analysis  of  a  dolphin  click. 


%  Filename :  dolphinspec . m 

%Author:  R.W.  Toland 

%Date  Last  Modified:  14  April  1998 

%Purpose:  This  program  calculates  the  energy  flux  spectral  density 
%  level  in  a  dolphin  click  and  plots  click+noise  and  noise 

%  prior  to  click  using  MATLAB  Normalization  (Parseval's 

%  Relation)  and  shows  that  the  time  plot  energy  equals  the 

%  frequency  plot  energy. 


clear  all 

%open  noise  data  file 

fidn  =  fopen ( Td: \toland\data\noiseprior 1 ) ; 

%read  data  file  that  is  in  binary  short  integer  format; 
sign  =  fread(fidn, inf , ’short* ) ; 

%normalize  the  signal 
sign  =  sign(l:128); 
signl  =  sign  -  mean(sign); 

Nndat  =  length (signl) ; 

%open  click  data  file 

fide  =  fopen  ( *  d: \toland\data\clickfc * )  ; 

%read  data  file  that  is  in  binary  short  integer  format; 
sige  =  f read (fide, inf, ’short* ) ; 

Inormalize  the  signal 

sige  =  sigc(l:128); 

sigcl  =  sigc(l:128)  -  mean(sigc); 

Ncdat  =  length (sigcl) ; 

% Sampling  Frequency 
Fs  =  1250000; 

%period  of  signal 
time  =  (l:Ncdat) /Fs; 

%convert  time  to  microseconds 
timel  =  time . *1000000; 

%plot  click 

figure (1) 

plot (timel, sigcl) 

title ( ’Bottlenosed  Dolphin  Echolocation  Click*) 
xlabel ( ’ Time  in  Microseconds  * ) 
ylabel ( *  Relative  Amplitude  * ) 

%fft  points 
Nfft  -  512; 

%calculate  fft  of  noise 
DFT  =  fft (sigcl, Nfft) ; 

Iconvert  data  points  to  frequency 
w  =  {0:{{Nfft/2)-l))/(Nfft/2)*(Fs/2); 

^convert  frequency  to  kHz 
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wl  =  w. /1000; 


%PSD  of  noise  prior  click 

[Pxxn,fn]  =  psd (signl, Nf ft, Fs, boxcar (Nfft)  ,  0)  ; 

%Multiply  by  2*Nfft/F s  to  get  true  PSD 
Pxxnl  =  Pxxn. * (2*Nfft/Fs) ; 

%Convert  to  dB  scale 

Pxxndb  -  10*logl0 (abs (Pxxnl) ) ; 

%PSD  of  click  +  Noise 

[Pxxc,fc]  =  psd ( sigcl, Nfft , Fs, boxcar (Nfft)  ,  0)  ; 

%Multiply  by  2*N/Fs  to  get  true  PSD 
Pxxcl  =  Pxxc. * (2*Nfft/Fs) ; 

%Convert  to  dB  scale 
Pxxcdb  =  10*logl0 (Pxxcl) ; 

N  =  length (Pxxcdb) ; 

%Change  units  of  Energy  Spectral  Density  PlotfdB  ref  lpicoJoule/mA2-Hz 
at  lm) 

fl  -  hydrocal (fn) ; 

%correction  for  hydrophone  (-dB  sensitivity  from  chart) 

Clickdb  =  Pxxcdb-fl; 

Noisedb  =  Pxxndb-fl; 

%correction  for  gain  (-201og (gain) ) ,  and  range  (+201og(d)),  correction 
for  energy  {-101og(pc)) 

Pxxcdb_new  =  Clickdb  -  14  +  20 . *logl0 ( 1 . 38 ) -10 . *logl0 ( 1 . 5e6) ; 
Pxxndb_new  =  Noisedb  -  14  +  20 . *logl0  (1 . 38 ) -10 . *logl0 (1 . 5e6) ; 

%Plot  both  (Click  +  Noise)  and  (Noise  prior  click)  on  same  plot 
figure (2) 

plot (fn (1:245) . /1000, Pxxcdbjiew (1:245) , T - 

1 z  fc ( 1 : 245)  ./1000, Pxxndb_new (1:245)/1-.1) 

title ( f Energy  Flux  Spectral  Density  Level  vs  Frequency’) 

xlabel ( T  Frequency ( kHz )  ? ) 

ylabel {' Energy  Flux  Spectral  .Density  Level  dB  re  [pJ/(mA2  Hz)]f) 
legend { 1 Click+Noise ? , ’Noise1 ) 

%Energy  in  frequency  domain  =  Energy  in  time  domain 
energy_from_time__click  =  sum  (sigcl .  A2 ) 
energy_from_freq__click  =  2*sum(Pxxc) 

energy_from_time__noise  =  sum  (signl .  A2 ) 
energy_from_freq_noise  =  2* sum (Pxxn) 
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%Filename :  hydrocal . m 
%Author:  R.W.  Toland 

%Date  Last  Modified:  14  April  1998 

%Purpose:  This  function  plots  the  ARL  430  hydrophone  sensitivity  curve 

%  out  to  the  first  600  Hz  to  apply  to  dolphinspec .m  for 

%  converting  PSD  plot  to  pressure  units.  The  curve  is 

%  approximated  by  polyfit  function. 

function  fl  =  hydrocal (fn) 

N  =  600; 

inta  =  round (225/600*N) ; 
intb  =  inta  +  round (75/600*N) ; 
intc  =  intb  +  round ( 80/600*N) ; 
intd  =  intc  +  round ( 70/600*N) ; 
inte  =  intd  +  round (70/600*N) ; 
intf  =  inte  +  round (30/600*N) ; 
intg  =  intf  +  round { 40/600*N) ; 
yl  =  -194 . *ones ( [inta, 1] ) ; 
for  j  =  1 : (intb-inta) 

y2(j)  =  yl (length (yl) )-j* (0.04)  ; 

end 

clear  j.; 

for  j  =  1: (intc-intb) 

y  3  ( j )  =  y2 ( length (y2) )-j*(0.013)  ; 

end 

clear  j; 

y4  =  y3 (length (y3) ) .*ones ([ (intd-intc) , 1] ) ; 
for  j  =  1: (inte-intd) 

y5(j)  =  y4 (length (y4 )) -j* (0. 0536)  ; 

end 

clear  j ; 

for  j  =  1: (intf-inte) 

y 6 ( j )  =  y5 (length (y5) )-j* (0.143) ; 

end 

clear  j ; 

for  j  =  1: (intg-intf) 

y 7  ( j )  =  y6 (length (y6) ) -j * ( . 03) ; 

end 

sum_y  =  [yl,,y2/y3ry4T/y5,y6,y7]  ; 
i  =  linspace (1, 600, length (sum_y) ) ; 
p  =  polyfit (i, sum_y, 9) ; 
f  =  polyval (p, (1: length (sum_y) )) ; 

fl  =  polyval (p, (fn. /1000) ) ; 
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APPENDIX  B.  CORRELATION  COEFFICIENT  ANALYSIS 


This  appendix  contains  the  MATLAB  program  used  to  conduct  correlation 
coefficients  and  comparison  analysis  of  five  consecutive  dolphin  clicks  from  a  click 
train. 


%Filename:  correlate. m 

%Author :  R.W.  Toland 

%Date  Last  Modified:  29  January  1998 

%Purpose:  This  program  loads  five  dolphin  clicks  containing  128  pts 
%  each,  and  correlates  each  click  to  one  of  the  other  five, 

%  producing  a  matrix  of  correlation  coefficients.  It  then 

%  plots  the  comparison  of  overlapping  signals. 


clear  all 
%open  data  file 

fida  =  fopen (’ d: \toland\data\clickl ')  ; 

%read  data  file  that  is  in  binary  short  integer  format; 
clicka  =  fread (fida, inf , ’short 1 ) ; 

%normalize  the  signal 

sigl  =  clicka  -  mean (clicka) ; 

N  =  length (sigl) ; 

%open  data  file 

fidb  =  fopen ( ’d: \toland\data\click2 1  )  ; 

%read  data  file  that  is  in  binary  short  integer  format; 
clickb  =  fread(fidb, inf,  'short 1 ) ; 

%normalize  the  signal 

sig2  =  clickb  -  mean (clickb) ; 

Nn  =  length(sig2) ; 

%open  data  file 

fide  =  fopen ( ' d: \toland\data\click3 1 )  ; 

%read  data  file  that  is  in  binary  short  integer  format; 
clickc  =  fread (fide, inf ,' short ')  ; 

■  %normalize  the  signal 
sig3  =  clickc  -  mean (clickc) ; 

Nc  =  length (sig3) ; 

%open  data  file 

fidd  =  fopen ( ?d: \toland\data\click4 1 )  ; 

%read  data  file  that  is  in  binary  short  integer  format; 
clickd  =  fread (fidd, inf, 'short1 ) ; 

%normalize  the  signal 

sig4  =  clickd  -  mean ( clickd) ; 

Nc  =  length (sig4) ; 

%open  data  file 

fide  =  fopen ( Td: \toland\data\click5 ’)  ; 

%read  data  file  that  is  in  binary  short  integer  format; 
clicke  =  fread (fide, inf , ’short'); 
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%normalize  the  signal 

sig5  =  clicke  -  mean (clicke ) ; 

Nc  =  length (sig5 ) ; 

% Sampling  Frequency 
Fs  =  1250000; 

%period  of  signal 
time  =  (1:N)/Fs; 

%convert  time  to  microseconds 
timel  =  time . *1000000; 

max  corr  =  [ ] ; 


for  i  =  1:5 

filel  =  [ ’ sig’ , num2str (i) ] ; 
for  j  =1:5 

file2  =  [ ' sig* , num2str ( j ) ] ; 

corr  =  xcorr (eval (filel) , eval (file2) , 1 coeff* ) ; 
[mx,indx]  =  max (abs (corr )) ; 
max__corr  (i,  j  )  =  mx; 
index (i,j)  =  indx; 
end 
end 


Correlatio.n_Coef  f  icient_Matrix  =  max_corr 
Index  Matrix  =  index 


for  k  =  1:129 
sigl_new (k+1, 
sig3_new (k+1, 
sig4_new (k+1, 
sig5__new  (k+1, 

end 


1)  =  sigl (k,  1 ) ; 
1)  =  sig3 ( k, 1 ) ; 
1)  =  sig4 ( k, 1 ) ; 
1)  =  sig5 ( k, 1) ; 


%plot  signal 
figure (1) 

plot  (timel, sigl_new (1:129),  timel,  sig2,  timel, sig3_new (1 : 129) , timel, 
sig4_new (1:129), timel, sig5_new (1 : 129) ) 
title ( ’Comparison  of  Five  Consecutive  Dolphin  Clicks1) 
xlabel(’Time  in  Microseconds’) 
ylabel ( 'Relative  Amplitude*) 


50 


APPENDIX  C.  SOUND  ABSORPTIVE  SCREEN  (SOAB)  ANALYSIS 


This  appendix  contains  the  MATLAB  programs  used  to  conduct  the  insertion 
loss,  reflection  loss  and  theory  analysis  of  SOAB. 


%Filename:  insertionloss .m 

%Author:  R.W.  Toland 

%Date  Last  Modified:  22  May  1998 

%Purpose:  This  program  calculates  the  insertion  loss  for  three 

%  different  sizes  of  SOAB  material. 

clear  all 

^frequencies  of  measurements 
freq  =  [130  240  250  270  290  340] ; 

%no  tile  hydrophone  receive  voltage 
Vo  =  [5.8  8.0  10.4  11.6  9.2  4.4]; 

Vol  =  [10.8  11.2  12.4  12.8  10.4  5.6]; 

Vo2  =  [1.5  7.2  9.4  14.0  14.4  6.0]; 

%with  tile  hydrophone  receive  voltage 
VI  =  [.7  1.5  1.75  2.1  2.1  1.45]; 

V2  =  [2.0  3.6  4.4  5.0  4.2  2.2]  ; 

V3  =  [.08  .5  .7  1.2  1.4  .7] ; 

V4  =  [.8  3.1  3.8  5.3  5.8  2.3]; 

%calculate  insertion  loss 
TL1  =  20 . *logl0 (Vo. /VI) ; 

TL2  =  20 . *logl0 (Vol . /V2 ) ; 

TL3  =  20.*logl0(Vo2./V3) ; 

TL4  =  20.*logl0(Vo2./V4) ; 

%plot  figure 

plot  (freq,  TL1,  freq,  TL1,  'o',  freq,  TL2 ,  '  :  ',freq,TL2,  ’  o ' ,  freq,  TL3,  ' 

freq,  TL3, 'o' , freq, TL4, freq, TL4 , 'o' ) 

title ('Tile  1  -  thin,  vacuum,  no  bubbles') 

xlabel ( ' frequency ( kHz ) ' ) 

ylabel (' Transmission  Loss (dB) ' ) 


%Filename:  insertionangle .m 

% Author:  R.W.  Toland 

%Date  Last  Modified:  22  May  1998 

%Purpose:  This  program  calculates  the  insertion  loss  as  a  function 

%  of  angle  off  axis. 

clear  all 

%frequency  of  measurements 

freq  =  [60  120  130  160  260  280  290  300  310  330  370  400] ; 

%hydrophone  recieve  voltage 

Vo  =  [2.45  4.1  3.7  1.7  2.85  3.6  3.7  3.3  2.9  2.55  1.9  1.5]; 

VI  =  [2.25  3.3  2.95  1.16  1.36  1.65  1.6  1.48  1.32  1.02  .66  .58]; 

V2  =  [2.25  3.2  2.9  1.25  1.36  1.65  1.6  1.4  1.28  1.02  .64  .58]; 

V3  =  [2.2  3.2  2.9  1.25  1.3  1.55  1.5  1.38  1.22  1.22  .64  .56]; 

V4  =  [2.2  3.15  2.85  1.2  1.25  1.5  1.45  1.35  1.2  .92  .6  .52]; 
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%calculate  insertion  loss 
TL1  =  20.*logl0(Vo./Vl) ; 

TL2  -  20 . *loglO (Vo . /V2 ) ; 

TL3  =  20  .  *logl0  (Vo.'/V3)  ; 

TL4  -  20.*logl0 (Vo./V4) ; 

%plot  figure 

plot (freq, TL1, freq, TLl,  ’ x ’ , f req, TL2 ,  1 : 1 , freq,TL3,  ’  .- 
1 ,  freq, TL4 , freq, TL4 ) 

title ( ’ SOAB  Insertion  Loss  As  A  Function  Of  Angle  Off  Axis’) 
xlabel ( ’ Frequency ( kHz )  ’  ) 
ylabel ( 1  Insertion  Loss(dB)’) 

legend ( ’SOAB’ ,  * ’ ,  ’SOAB-5deg’ ,  ’SOAB-lOdeg’ ,  ’ SOAB-20deg’ ) 


%Filename :  iltheory .m 
%Author :  R.W.  Toland 

%Date  Last  Modified:  22  May  1998 

% Purpose :  This  program  calculates  the  insertion  loss  for  1/2"  SOAB 
%  filter  and  compares  it  to  theory. 

clear  all 

freq  =  [60  130  260  300  330  400].*le3; 

%hydrophone  receive  voltage 
Vo  =  [4.2  6.0  3.2  4.1  3.7  2.0]; 

V2  =  [3.8  4.5  1.4  1.6  1.28  0.72]; 

TL2  =  20 . *logl0 (Vo . /V2 ) ; 

%plot  measured  values 

plot (freq. /1000, TL2,  ’ - ’ , f req. /1000,  TL2 ,  ?x’) 
title (’ Insertion  Loss  Of  1/2"  SOAB  Tile’) 
xlabel ( ’ frequency ( kHz ) ’ ) 
ylabel { ’ Insertion  Loss (dB) ’ ) 

%grid 
hold  on 

p  =  polyfit (freq, TL2, 3) ; 

freq_theory  =  [60  130  260  300  330  400]-*le3; 

TL_theory  =  polyval (p, freq_theory)  ; 

%Insertion  Loss  Theory 
w  =  2*pi^freq_theory; 


pi  =  1000; 
cl  =  1480; 


%density  fresh  water 
%speed  of  fresh  water 


p2  =  1662; 

c2  -  933.33; 
d  =  .5/39.4; 
k  =  w. /c2; 


%density  SOAB 

%longitudinal  sound  speed  SOAB 
%thickness  in  inches  converted  to  meters 


% longitudinal  attenuation 

alpha_l  —  TL_theory . * 100 . /8 . 7 ;  %conversion  of  dB/cm  to  nepers/m 

%longitudinal  loss  parameter 
r  =  alpha_l*c2 . /w; 
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Z1  =  (pl*cl); 

Z2  =  (p2*c2 ) ; 

m  =  Z2./Z1; 

num  =  ( (m+l-j . *r) . A2 . *exp ( j*2 . *k*d) ) - { (m-l+j . *r) . A2 . *exp (- 
2 . *alpha_l*d) ) ; 

denom  =  4*m* ( 1-j . *r) . *exp ( ( j . *k*d) - (alpha_l*d) ) ; 

IL  =  20*logl0 (abs (num. / denom) ) ; 

%plot  theory 

plot (freq_theory. /1000, IL,  *  - . 1 ) 
gtext ( ’Measured  (-) 1 ) 
gtext ( ’Theory  (-. ) ’ ) 


%Filename:  rltheory .m 
%Author:  R.W.  Toland 

%Date  Last  Modified:  22  May  1998 

%Purpose:  This  program  calculates  the  reflection  loss  for  1/2”  SOAB 

%  filter  and  compares  it  to  theory.  . 


clear  all 

freq  -  [80  130  260.  300  330  400]; 

%SOAB  III  (1/2") 

Vi  =  [.72  1.6  4.2  8.2  9.0  1.7]; 

Vr  =  [  .034  .08  .42  .76  .9  .2] ; 

RL  =  20 . *logl0 (Vi . /Vr) ; 
figure (1) 

%plot  measured 

plot ( freq, RL, freq, RL, ’ x ’ ) 

title ( ’Reflection  Loss  (1/2”  SOAB  Tile)’) 

xlabel ( ’ Frequency  ( kHz ) ’ ) 

ylabel (’ Ref lection  Loss  (dB)’) 

axis ( [50  400  0  32]  ) 

set (gca, ’ ydir ’ , ’ reverse ’ ) 

hold  on 

freq_theory  =  [ 60 : 10 : 400] . *le3; 

%Reflection  Loss  Theory 
w  =  2*pi*freq_theory; 
pi  =  1000;  %density  fresh  water 

cl  =  1480;  %speed  of  fresh  water 

p2  =  993.33;  %density  SOAB 

c2  =  1662;  %longitudinal  sound  speed  SOAB 
d  -  .5/39.4;  %thickness  in  inches  converted  to  meters 
k  =  w. /c2 ; 

%longitudinal  attenuation  (freq  in  kHz) 
a  =  0.20;  %Np/m  kHZ 

b  =  12e-5;  %Np/m  kHz 

alpha_l  =  (a. *freq_theory . /1000 )  +  (b. * ( freq_theory . /1000 )  . A2)  ; 
%longitudinal  loss  parameter 
r  =  alpha_l*c2 . /w; 
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Z1  =  (pl*cl) ; 

Z2  =  (p2*c2 ) ; 
m  =  Z2./Z1; 

numl  =  (  (m+l-j . *r) . / (m-l+j . *r) ) . *exp ( ( j . *k*d)  +  (alpha_l*d) ) ; 
mm2  =  (  (m-l+j  .  *r)  .  /  (m+l-j  .  *r)  )  .  *exp  (  (-j  .  *k*d)  -  (alpha_l*d)  )  ; 
denom  =  exp  ( ( j . *k*d)  +  (alpha_l*d) ) -exp ( (-j . *k*d) - (alpha_l*d) ) ; 
ER  =  20*logl0 (abs ( (numl-num2) . /denom) ) ; 

%plot  theory 

plot (freq_theory . /1000, ER,  ' - . 1  ) 
set (gca, ' ydir 1 , f  reverse T ) 
gtext ( ’Measured  (solid) ’ ) 
gtext ( ’Theory  (dashed)  ’ ) 
gtext ( ’ a  =  0.20’) 
gtext (*b  =  12*le-5’) 
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APPENDIX  D.  WPLOT  PROGRAM  SAMPLE 


This  appendix  contains  an  example  of  the  Wplot  software  used  to  visually  display 
raw  data  echolocation  recordings.  This  software  allowed  large  data  sets  of  interest  to  be 
cut  into  manageable  files  which  could  then  be  analyzed  in  MATLAB.  [Ref.  20] 
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